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We discuss the production of ultracold molecules in their electronic ground state by photoassocia- 
tion employing electronically excited states with ion-pair character and strong spin-orbit interaction. 
A short photoassociation laser pulse drives a non-resonant three-photon transition for alkali atoms 
colliding in their lowest triplet state. The excited state wave packet is transferred to the ground 
electronic state by a second laser pulse, driving a resonant two-photon transition. After analyzing 
the transition matrix elements governing the stabilization step, we discuss the efficiency of popula- 
tion transfer using transform-limited and linearly chirped laser pulses. Finally, we employ optimal 
control theory to find the most efficient stabilization pathways. We find that the stabilization ef- 
ficiency can be increased by one and two orders of magnitude for linearly chirped and optimally 
shaped laser pulses, respectively. 

PACS numbers: 82.53.Kp,33.80.-b,31.50.-x,33.90.+h 



I. INTRODUCTION 

Photoassociation, forming molecules from ultracold 
atoms using laser light [1], is a prime candidate for coher- 
ent control which utilizes the wave nature of matter in 
order to steer a process such as bond formation toward a 
desired target [2, 3]. At very low temperature, the deli- 
cate build-up of constructive and destructive interference 
between different quantum pathways is not hampered by 
thermal averaging. The basic tool for coherent control 
are short laser pulses that can be shaped in their ampli- 
tude, phase and polarization. They can drive both adia- 
batic and non-adiabatic photoassociation dynamics [4]. 

A particular feature of photoassociation at very low 
temperatures is the excitation of an atom pair at fairly 
large interatomic separations [1]. This results from a 
compromise between the atom pair density in the elec- 
tronic ground state, highest at large interatomic separa- 
tions, and population of excited state bound levels with 
reasonable binding energies, that increase with decreas- 
ing interatomic separations. Therefore, the frec-to-bound 
transition matrix elements are largest for photoassocia- 
tion at separations of 50 ao to 150 ao with corresponding 
detunings of less than about 20 cm^ 1 . Although these 
matrix elements are optimally chosen, they are several 
orders of magnitude smaller than those for the excita- 
tion of atoms. This poses a problem for photoassoci- 
ation with short laser pulses which inherently have a 
large bandwidth. As soon as the wings of the pulse spec- 
trum overlap with the atomic resonance, atoms instead 
of bound levels are excited [5], and subsequent sponta- 
neous emission depletes the trapped sample [6, 7]. In 
photoassociation experiments using broadband femtosec- 
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ond laser pulses, the pulse spectrum therefore needed to 
be cut to suppress excitation of the atomic resonance [8] . 
The sharp spectral cut yields long wings of the tempo- 
ral pulse profile, and the ensuing photoassociation dy- 
namics were dominated by transient Rabi oscillations of 
extremely weakly bound molecules caused by the long 
tail of the pulse [8-11]. While it was gratifying to see 
that femtosecond photoassociation is feasible [8], larger 
binding energies and vibrational instead of electronic dy- 
namics are required to produce stable molecules in their 
electronic ground state [4]. 

Femtosecond photoassociation at very low tempera- 
ture requires driving a narrow-band transition with a 
broad-band laser. This can be achieved by employing 
multi-photon rather than one-photon transitions [12]. 
The high peak powers of femtosecond laser pulses eas- 
ily allow for driving multi-photon transitions, and multi- 
photon control schemes have been demonstrated for both 
weak [12, 13] and strong laser pulses [14-16]. In the weak- 
field regime, perturbation theory shows that optical in- 
terference of two or more photons can be used to com- 
pletely suppress absorption [12, 13]. For intermediate 
intensities, higher order perturbation theory can be em- 
ployed to obtain rational pulse shapes [17-21] that allow 
to control the absorption. In the strong-field regime, dy- 
namic Stark shifts drive the transition out of resonance. 
This can be countered by a linear chirp of the pulse which 
compensates the phase accumulated due to the Stark 
shift. Adjusting the amplitude of the pulse to guarantee 
a 7r or 2-7T pulse controls the absorption [14-16]. These 
control schemes can be applied to femtosecond photoas- 
sociation in order to suppress the excitation of atoms [22]. 
Multi-photon transitions can also be utile for femtosec- 
ond photoassociation at high temperature [23, 24]. There 
the main advantage derives from the larger flexibility 
in transition energies, obtained when combining two or 
more photons, and the new selection rules. The disadvan- 
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tage of high temperatures is the low initial coherence, or 
quantum purity, of the thermal ensemble of atoms. How- 
ever, femtosecond photoassociation can generate rovibra- 
tional coherence by Franck-Condon filtering [23]. The 
long-standing goal of coherently controlling bond forma- 
tion [25] has thus become within reach also for high tem- 
peratures. 

Besides the possibility of driving a narrow-band tran- 
sition, multi-photon femtosecond photoassociation also 
allows for accessing highly excited electronic states that 
may have significant ion-pair character. Such states are 
expected to be well suited for the formation of stable 
molecules in their electronic ground state due to the 
peculiar shapes of the potential energy curves obtained 
when an ion-pair state crosses covalent ones [26, 27]. 

These effects become most significant for heavy atoms 
with strong spin-orbit interaction. The coupling of two 
(or more) electronic states leads to strong mixing of the 
rovibrational levels provided the coupling becomes reso- 
nant [28] . The wavefunctions of such strongly mixed lev- 
els display peaks at all the four classical turning points 
providing large transition matrix elements for both pho- 
toassociation and subsequent stabilization to the elec- 
tronic ground state [29]. For homonuclcar diatomics, 
usually several neighbouring vibrational levels are af- 
fected by the resonant coupling [30, 31], making them 
particularly suitable for short pulse photoassociation and 
stabilization since a laser pulse addresses a wave packet, 
not just a single level. In the case of heteronuclear 
molecules, the resonantly perturbed levels are typically 
isolated within the vibrational spectrum. However, the 
peaks at the inner turning points are so large that sta- 
bilization into deeply bound levels of the ground state 
well [32, 33] all the way down to v" = for SrYb [34] 
becomes feasible in a single step. Strong spin-orbit in- 
teraction furthermore allows for singlet-triplet conver- 
sion [35, 36]. 

Here, we combine all these features in a study of short- 
pulse multi-photon photoassociation into highly excited 
states with significant ion-pair character and strong spin- 
orbit interaction. Our envisioned scheme for off-resonant 
three-photon photoassociation and subsequent resonant 
two-photon stabilization is displayed in Fig. 1. The po- 
tential energy curves shown in Fig. 1 as well as the spin- 
orbit couplings and transition matrix elements were cal- 
culated with state of the art ah initio methods. The inter- 
action of the atom pair with the laser pulses is modelled 
non-perturbatively. A photoassociation pulse excites a 
pair of rubidium atoms, colliding in their lowest triplet 
state, into the manifold of electronically excited states 
just below the 5s + 6p 1 / 2 threshold. The shape of the po- 
tential energy curves close to the Franck-Condon radius 
for photoassociation is dominated by the 1 /R behavior of 
the ion-pair state. Note that excitation into the ion-pair 
state is dipole-forbiddcn. The three-photon photoasso- 
ciation indicated in Fig. 1 is facilitated by the covalent 
states that are coupled to the ion-pair state. Upon its cre- 
ation by the photoassociation pulse, the wave packet rolls 
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FIG. 1: (Color online) Proposed scheme for multi-photon pho- 
toassociation and subsequent stabilization producing ultra- 
cold Rb2 molecules their electronic ground state. 



down the 1/R slope, being subject to spin-orbit oscilla- 
tions between different electronic states of the manifold 
at the same time. Once the wave packet reaches short 
interatomic separations, a stabilization pulse catches it, 
transferring it to the electronic ground state via a reso- 
nant transition into the intermediate A 1 Yi+-b 3 U u states. 
The strong spin-orbit coupling in both the highly ex- 
cited and intermediate state affords a conversion of the 
molecules to purely singlet character in the final state. 
Our emphasis in the present study is on the intricate 
dynamics of the photoassociated wave packet in the elec- 
tronic manifold below the 5s + 6pi/2 asymptote and the 
stabilization to the electronic ground state. We employ 
rationally shaped laser pulses as well as optimal control 
theory to determine the maximum stabilization efficien- 
cies, taking standard constraints of pulse shaping exper- 
iments into account. The application of optimal control 
theory to the photoassociation step requires a theoret- 
ical description that accounts for the initial incoherent 
ensemble of atoms and will be presented elsewhere. 

The plan of our paper is as follows. The theoretical 
framework is presented in Sec. II, describing the model 
in Sec. II A, details of the electronic structure calcula- 
tions and optimal control theory in Sec. II B and II C, 
respectively. The creation of the photoassociated wave 
packet, representing the initial state for the stabilization 
step, and its dynamics in the coupled manifold of elec- 
tronic states below the 5s + 6pi/2 asymptote is studied in 
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Sec. III. The stabilization efficiency of transform-limited 
and linearly chirped pulses is investigated in Sec. IV. 
Optimally shaped pulses driving the excited state wave 
packet into the electronic ground state and their effi- 
ciency are discussed in Sec. V. We conclude in Sec. VI. 

II. THEORETICAL FRAMEWORK 

A. Hamiltonian 

We consider a pair of 85 Rb atoms, held at a temper- 
ature of 100 fiK, typical for magneto-optical traps, col- 
liding in the a 3 E+ lowest triplet state. The formation 
of molecules by photoassociation and stabilization of the 



excited state molecules to the eletronic ground state are 
treated separately. First, a photoassociation laser pulse 
drives a three-photon transition, red detuned with re- 
spect to the 2 S(5s) + 2 Pi/2(6p) asymptote, creating a 
molecular wave packet in the manifold of the (5) 1 S+, 
(6)^, (7yj:j, (3) 3 n 5 , and (4) 3 n 5 electronically ex- 
cited states that partially have an ion-pair character, cf. 
Fig. 2. The states in this manifold are coupled by the 
spin-orbit interaction and by the non-adiabatic radial 
coupling matrix elements. In our calcuations, only the 
(6) 1 !];]", and (3) 3 n g components of the photoas- 
sociated wave packet turned out to be significant. Ne- 
glecting the (7) 1 !];*" and (4) 3 II S states, the Hamiltonian 
describing the three-photon photoassociation reads 
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in the (three-photon) rotating wave approximation. In 

„ 2S+1 |A| 

Eq. (1), H denotes the Hamiltonian for the nuclear 

motion in the 25+1 |A| electronic state, 

H 2S+1|A| =t + V 2S+1 ^(R) + 4 S+1 ^(t,R) + A UL , (2) 

with the kinetic energy operator given by T = ^ , /i 

the reduced mass and V" S+1 l A l(i?) the potential energy 
curve. The three-photon detuning, , is taken with re- 
spect to the atomic 2 S(5s) — > 2 P\/2{'op) three-photon 
transition. For strong photoassociation laser pulses, the 

2S+1 1 Al 

dynamic Stark shift, uj s (t,R), will become signifi- 
cant. It arises from the interaction of the 2S+1 |A| state 
with the intermediate off-resonant states and is given 
by the effective dynamic electric dipole polarizability, 

a e ff(u>L,R), 

4 S+1|A ' (t, R) = -i|e(t)| 2 a e// K, R) , (3) 

where e(t) = |e(i)|e 11 ^*) describes the electric field of 
the laser pulse in the rotating frame with the envelope 



\e(t)\ and (f>(t) denoting the relative phase, taken with 
respect to the central frequency's phase. x' 3 '(u;£,i?) is 
the three-photon electric dipole transition moment, £,i(R) 
(i = 3, 5) are the spin-orbit couplings and A(R) is the 
non-adiabatic radial coupling matrix element between 
the (5)^+ and (G) 1 !^ states. 



In a second step, the initial wave packet created by the 
three-photon photoassociation, is deexcited to the A^E^ 
ground electronic state via a resonant two-photon elec- 
tric dipole transition. The intermediate states for the 
two-photon transition are the A 1 ^ and 6 3 n„ states, 
correlating to the 2 S(5s) + 2 P(5p) asymptote, that are 
also strongly coupled by spin-orbit interaction. Electric 
dipole transitions are allowed between all components 
of the initial wave packet and the intermediate states, 
whereas the A' 1 E+ ground electronic state is only con- 
nected to the ^4 1 Sj state by a strong electric dipole tran- 
sition. The Hamiltonian describing the stabilization of 
the photoassociated wave packet to the electronic ground 
state via a resonant two-photon transition reads 
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in a (one-photon) rotating wave approximation. The 
Hamiltonian for field-free nuclear motion in the 25+1 |A| 

- 2S+1 |A| 

electronic state, H , is now given by 

h 2S+1|a| = t + y 2S+1 i A i(i?) + A™p (5) 

with the dctvmings A^f = W2p( 5p ) — ljl and = 
W2p( 6p ) — 2ujl for the states dissociating into the 2 5(5s) + 
2 P(5p) and 2 S(5s) + 2 P(6p) asymptotes, respectively. 
The electric transition dipole moments are denoted by 

d 1 (R) = (X 1 E+|a| J 4 1 S+) , 
d 2 (R) = {A l E+\d\(5)^), 

d 3 (R) = (6 3 n u |d|(3) 3 n g ), 

d 4 (R) = (A 1 E+|a|(6) 1 S+) J 
and the spin-orbit coupling matrix elements read 

UR) = (& 3 n u |H5=°Vn M ), 

UR) = ((5) 1 E+|H^|(3) 3 n g ), 

UR) = ((3) 3 n ff |H^|( 3 ) 3 n 9 ), 
= ((e^s+lH^ |(3) 3 n ff ), 

where Hso denotes the spin-orbit Hamiltonian in the 
Breit-Pauli approximation. For large interatomic sepa- 
rations R, the transition dipole moments and spin-orbit 
coupling approach their atomic values, 

di(i?->oo) = \/2( 2 S'(5s)|d| 2 P(5p)) , 
d 2 (R^ oo) = 0, 
d 3 (R^ oo) = 0, 

d 4 (R ->oo) = V2( 2 5(5s)|d| 2 F(5p)) , 

and 

oo) = £i(fl->oo)/\/2 = A FS ( 2 P(6p))/3, 
£ 4 (i?^oo) = &(R^ oo)/V2 = A fs ( 2 P(5jj))/3, 
C 5 (i?^^) = 0, 
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FIG. 2: (Color online) The {5) 1 E+, (6)^+, and (3) 3 n s elec- 
tronic states relevant for initial wave packet dynamics and 
ion-pair curve. 

where A FS { 2 P{5p)) = 237.6cm" 1 and A FS ( 2 P(6p)) = 
77.5 cm -1 are the atomic fine structure splittings. 

The Hamiltonians, Eqs. (1) and (4), are represented 
on a Fourier grid with an adaptive step size [37-39] using 
N = 2048 grid points. The time-dependent Schrodinger 
equation for the pump and dump Hamiltonians, 

i^|*(t)) = H|*(t)>, (6) 
is solved with the Chebyshev propagator [40]. 

B. Electronic structure 

State-of-the-art ab initio techniques have been applied 
to determine the electronic structure data of the rubid- 
ium molecule needed in our dynamical models of the 
pump and dump processes. The details of the calcu- 
lations and a thorough discussion of the accuracy by 
comparison with the most recent high-resolution spec- 
troscopic results will be reported elsewhere [41]. Briefly, 
all potential energy curves for the gcrade and unger- 
ade singlet and triplet states lying below 26000 cm -1 
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FIG. 3: (Color online) Electric transition dipole moments (a) 
and spin-orbit coupling matrix elements (b) essential for stud- 
ied dynamics. 



at the dissociation limit were calculated with the re- 
cently introduced Double Electron Attachment Interme- 
diate Hamiltonian Fock Space Coupled Cluster method 
restricted to single and double excitations (DEA-IH-FS- 
CCSD) [42, 43]. Starting with the closed-shell refer- 
ence state for the doubly ionized molecule Wa^ + that 
shows the correct dissociation at large distances R into 
closed-shell subsystems, Rb + +Rb+, and using the dou- 
ble electron attachment operators in the Fock space cou- 
pled cluster ansatz renders our method size-consistent 
at any interatomic distance R and guarantees the cor- 
rect largc-i? asymptotics. Thus, the DEA-IH-FS-CCSD 
approach overcomes the problem of the standard CCSD 
and equation of motion CCSD methods [44] with the 
proper dissociation into open-shell atoms. The poten- 
tial energy curves obtained from the ab initio calcula- 
tions were smoothly connected at intermediate distances 
to the asymptotic multipolc expansion [45]. The Cg co- 
efficient of the electronic ground state and the C3 coeffi- 
cient of the first excited state were fixed at their empirical 
values derived from high-resolution spectroscopic experi- 
ments [46, 47], while the remaining coefficients were taken 
from Ref. [48]. 

Electric transition dipole moments, radial non- 
adiabatic coupling, and spin-orbit coupling matrix ele- 
ments were obtained by the mult i- reference configura- 
tion interaction method (MRCI) , restricted to single and 
double excitations with a large active space. Scalar 
relativistic effects were included by using the small- 
core fully relativistic energy-consistent pscudopotcntial 
ECP28MDF [49] from the Stuttgart library. Thus, in the 
present study the Rb2 molecule was treated as a system 
of effectively 20 electrons. The [14sl4p7d6/] basis set was 
used in all calculations. This basis was obtained by de- 
contracting and augmenting the basis set of Ref. [49] by a 
set of additional functions, improving the accuracy of the 



atomic excitation energies of the rubidium atom with re- 
spect to the NIST database. With this basis set, employ- 
ing the DEA-IH-FS-CCSD method for the nonrelativistic 
energies and the MRCI approach for the fine structure, 
we reproduce the experimental excitation energies with 
a root mean square deviation (RMSD) of 21 cm" 1 . 

The DEA-IH-FS-CCSD calculations were done with 
the code based on the ACES II program system [50], 
while the MRCI calculations were performed with the 
MOLPRO code [51]. 

The accuracy of the computed potential energy curves 
is as excellent as that of the atomic excitation ener- 
gies [41]. Here, we only point out that the RMSD of 
the well depths of the present potential energy curves 
for electronic states that were observed experimentally 
is 64 cm" 1 . In particular, the well depths of the ground 
singlet state and the lowest triplet state are reproduced 
within 2.7% and 3.7%, respectively. The relative er- 
rors for the excited states relevant for the present study, 
the A X T,^ and b 3 H u pair, are even smaller, 0.5% and 
0.3%, respectively. We expect a similar accuracy for the 
higher electronic states that have not yet been observed 
experimentally. The accuracy of the present results for 
the potential energy curves is much higher than that of 
Refs. [52, 53] and slightly better than in the recent study 
by Allouche and Aubert-Frecon [54] which did not con- 
sider the electronic states crucial for our photoassociation 
proposal. 



C. Optimal control theory 

Optimal control theory (OCT) can be used to calculate 
the shape of laser pulses that efficiently drive a desired 
transition. We will employ it here to determine the most 
efficient stabilization between an initial molecular wave 
packet and deeply bound levels in the ground electronic 
state. In principle, this problem is compeletely control- 
lable such that perfect population transfer can be real- 
ized. However, contraints such as limited pulse duration, 
spectral bandwidth and pulse intensity will compromise 
the stabilization process, reducing the transfer efficiency. 

The control problem is defined by minimization of the 
funcional 



J = J T 



9W)]dt, 



(7) 



where the first term denotes the final-time T target and 
the second one intermediate-time costs. The final-time 
target, Jt, for a state-to-state transition from an initial 
state \^i n ) to a target state, here a vibrational level of 
the electronic ground state, v" , is written as 



J™ = 1-\($ V „\U(T,0; c)|*< 



(8) 



U(T, 0; e)!^^) represents the formal solution of the time- 
dependent Schrodinger equation, i.e., J ss corresponds 
the overlap of the initial state, propagated to the final 
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time T under the action of the laser field e(t), with the 
target state. Instead of optimizing the transition to a 
single vibrational level, we can also choose the transition 
into a manifold of states for the target. This is expressed 
by the functional 



J£ m = 1 - 



K#„«|U(T, 0;e)\% 



(9) 



where any vibrational level of the electronic ground state 
with quantum number between u^ in and u^ 1QX can be 
populated at final time. Once the optimum is reached, 
both functionals, J ss and J sm , take the value zero. 

The intermediate time cost, g[e(i)], can in general de- 
pend on both the state and the field. Here we restrict 
the dependence to one on the laser field only, 



9[<t)] 



A 

W) 



E 

a={9U,3m} 



(*) 



(t) 



(10) 



This choice minimizes the change in the pulse energy be- 
tween iterations k and fc+1 [55]. The shape function S(t). 
S(t) = sin 2 (nt/T), enforces a smooth switch on and off 
of the field and A is a weight. Note that the laser field is 
complex since we employ the rotating-wave approxima- 
tion. A non-zero phase indicates a relative phase with 
respect to the laser pulse peak center, or, in the spectral 
domainc, with respect to the phase of the central laser 
frequency. 



Using the linear variant of Krotov's method [56, 57], 
the update equation for the laser field at iteration step 
k + 1 can be derived, 
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where U(i, 0; e^" 1 " 1 ))^^) is the initial state forward 
propagated to time t under the action of the new field, 
e< fc+1 ), and (i(t, T; e( fc ))|^„«) denotes the target state(s) 
backward propagated to time t under the action of the 
old field, e( fe ). The derivative of the Hamiltonian with 
respect to the field yields a matrix having as its elements 
all the transition dipole moments di, cf. Eq. (4). Opti- 
mization of the functionals J ss or J sm requires repeated 
forward and backward propagations of the initial and tar- 
get states. 



III. EXCITED STATE WAVE PACKET 
REPRESENTING THE INITIAL STATE FOR 
STABILIZATION 

The initial wave packet for the stabilization step is cre- 
ated by the photoassociation pulse. The simplest pulse 
that can be employed for the three-photon photoassoci- 
ation is a transform-limited (TL) Gaussian pulse. The 
intensity of the laser pulse is chosen to be in the per- 
turbative week-field regime, where the composition of 
the photoassociated wave packet reflects the bandwidth 
of the laser pulse combined with the vibrationally av- 
eraged three-photon electric dipole transition moments 
between the initial scattering state and the excited state 



vibrational levels below 2 S(5s) + 2 P\/2{&p) dissociation 
limit. A pulse duration of 4 ps full width at half max- 
imum (FWHM) is considered, corresponding to a spec- 
tral bandwidth of 3.7 cm -1 . The pulse is red detuned by 
12 cm -1 from the 2 S(5s) — > 2 -f\/2(6p) atomic three- 
photon transition. In order to utilize broadband fem- 
tosecond laser pulses, more elaborate pulse shapes will be 
required that suppress the excitation of atoms [22] while 
possibly maximizing frce-to-bound transitions. However, 
the general features of the photoassociated wave packet 
are determined by the three-photon matrix elements. 
They are the largest close to the 2 S*(5s) + 2 P\/2{ < op) 
dissociation limit, corresponding to a photoassociation 
window at interatomic separations between 30 and 
50 ao. Transition moments to the vibrational levels de- 
tuned by more then 30 cm -1 from the 2 5'(5s) + 2 Pi/2(6p) 
threshold, corresponding to photoassociation windows at 
shorter interatomic separations, are significantly smaller. 
The most important contributions to the photoassoci- 
ated wave packet will therefore remain the same as in our 
simple example. It could turn out that, using coherent 
control, photoassociation into wave packets with binding 
energies larger than 30 cm -1 becomes feasible. In this 
case, the stabilization of the photoassociated molecules 
becomes easier, and their faster vibrational dynamics and 
larger Franck-Condon factors to the deeply bound A rl S+ 



7 



0.08 



0.06 



3. 0.04 



0.02 



0.00 




20 30 
R (bohr) 



50 



FIG. 4: (Color online) Wave packet obtained after photoas- 
sociation with a transform-limited pulse. The pulse duration 
is 4ps (FWHM) and a snapshot 2.5 ps after the photoasso- 
ciation pulse maximum is shown. At this time, the (6) 1 S ! J" 
component is insignificant and not visible on the scale of this 
figure. Inset: Decomposition of the wave packet onto the vi- 



brational levels of the coupled (5) 1 E ! j 
states. 



(6)^+, and (3) 3 n 9 



vibrational levels will only improve the predictions of the 
present study. 

The initial wave packet for the stabilization step is 
plotted in Fig. 4. The inset shows the decomposition 
of the wave packet onto the vibrational levels of the cou- 
pled (6) x i:+, and (3) 3 n 3 states. The binding 
energy of the wave packet amounts to 11.55 cm -1 . The 
snapshot shown in Fig. 4 is taken 2.5 ps after the maxi- 
mum of the photoassociation pulse, i.e., before the pulse 
is over. At this time, the Gaussian character of the wave 
packet is still apparant, while at later times the strongly 
anharmonic shape of the potential energy curves leads to 
strong wave packet dispersion. Note that the photoas- 
sociated wave packet shows truly mixed character with 
about 65% of its norm residing on the (3) 3 IT g triplet 
component and 35% on the (S) 1 !]^" singlet component. 
This is despite the fact that electric dipole transitions 
are allowed only between the atomic pair in the a 3 S+ 
lowest triplet state and the (3) 3 II fl triplet component of 
the coupled electronically excited manifold dissociating 
into 2 S(5s) + 2 -Pi/2(6p) and illustrates the strong spin- 
orbit coupling. The role of the spin-orbit coupling is 
further evidenced by the double peak structure of the 
wave packet with the short-range peak corresponding to 
the outer turning point of the (3) 3 FI ff potential and the 
long-range peak reflecting the outer turning point of the 
(5) 1 E+ potential. The importance of the strong resonant 
spin-orbit coupling for the stabilization of photoassoci- 
ated molecules has been discussed extensively [29-34, 58]. 
In the present study, we will not only use it for improved 
transition matrix elements to deeply bound ground state 
levels, but also to convert an atom pair colliding in the 
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FIG. 5: (Color online) a: Time evolution of the populations 
of the (5) 1 E^ _ , (6) 1 S^", and (3) 3 I1 9 components of the wave 
packet, b: Time evolution of the average bond length of the 
wave packet. 



triplet state to singlet molecules [35]. 

The initial wavcpacket propagates toward shorter in- 
teratomic separations under the influence of the excited 
state potentials. At large interatomic separations, the 
potential energy curve of the state displays a 

strong —1/R ion-pair character. The singlet-triplet oscil- 
lations are analyzed in Fig. 5(a) displaying the singlet and 
triplet components of the wave packet evolving after the 
photoassociation pulse in the manifold of electronically 
excited states. The population of the (3) 3 IT ff triplet com- 
ponent oscillates around 60%, whereas the population of 
the (6) 1 E J J" component, that was absent just after pho- 
toassociation, reaches a maximum of 27% at t = 12.7 ps 
after the peak of photoassociating pulse. A second maxi- 
mum of the (6) 1 E+ component is observed after a period 
of 20.1 ps and a third one after another 36.2 ps. The times 
at which the (6) 1 E^j" component reaches maximal values 
can be interpreted as moments when the wave packet 
arrives at its shortest distance and is reflected from the 
innermost repulsive short range wall. This observation is 
confirmed by calculating the average bond length of the 
wave packet, shown in Fig. 5(b), which allows to estimate 
the revival time of the present wave packet to be between 
20 ps and 30 ps. This estimate agrees with the range of 
revival times, defined by T rev (v) = i-g — xtP — _ 2E i [59], 



for the vibrational levels v that make up a wave packet 
with binding energies close to 12 cm -1 . 

The knowledge of the revival time of the wave packet 
is useful for the interpretation of the projections of the 
time-dependent wave packet onto vibrational levels of 
the ^ 1 E+ ground electronic state, shown in Fig. 6(a) 
for all ground state levels with binding energies up to 
1000 cm" 1 . These projections are largest when the wave 
packet is localized at its inner turning point, cf. Fig. 5(b). 
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FIG. 6: (Color online) a: Projection of the time-dependent 
wave packet, and its (5) 1 E+, (6)^+, and (3) 3 I1 9 compo- 
nents, onto all vibrational levels of the X 1 H'£ ground elec- 
tronic states with binding energies up to 1000 cm -1 . b: Pop- 
ulation of the X 1 ^ component of the wave function after 
stabilization with a TL pulse vs time delay between photoas- 
sociation and stabilization pulse for three different detunings, 
taken with the respect to the 2 -Pi/2(6p) — > 2 S(5s) two- 
photon transition. The stabilization pulse has a FHWM of 
100 fs and an integrated pulse energy of 25.4 nj correspond- 
ing to the weak field regime. The arrow indicates time 13.3 ps 
chosen as the the center of pulses. 



The time at which the projections show maxima corre- 
spond to optimal time delays between photoassociation 
and stabilization pulse. The times inbetween these max- 
ima are given by the revival time. However, the transi- 
tion probability does not only depend on the overlap of 
initial and final wave function, but also on the dipole mo- 
ments and the topology of the intermediate state surfaces 
and their coupling. This is illustrated by the difference 
between Fig. 6(a) and (b) with (b) showing the calcu- 
lated population on the X 1 ^ ground electronic state 
as a function of the time delay between photoassociation 
and stabilization pulse for three different detunings of the 
stabilization pulse. The -X -1 ^ population is obtained by 
solving the time-dependent Schrodinger equation for a 
weak TL stabilization pulse. It is often argued that a 
weak pulse projects the wave packet onto the final state. 
Such a simple picture is not valid here where manifolds of 
coupled states rather than single isolated electronic states 
are present. The positions of the maxima and minima of 
the final X 1 Yj+ state population correspond to those of 
the total projection (black line in Fig. 6(a)), but the de- 
tails agree much better with the time-dependence of the 
projection of the (6) 1 S^J" component (dashed blue line in 
Fig. 6(a)). This suggests that the transition from the 
(6) 1 S^ _ component is the most important one in the sta- 
bilization process. 

Based on the time-dependence of the projection and 
the stabilization probability analyzed in Fig. 6, we choose 



to = 13.3 ps for the time delay, taken between the peak 
of the photoassociating pulse and the center of all pulses 
used in the following sections. 



IV. STABILIZATION TO THE ELECTRONIC 
GROUND STATE WITH TRANSFORM-LIMITED 
AND LINEARLY CHIRPED PULSES 

We first study transform-limited and linearly chirped 
stabilization pulses, in order to understand the role of 
the basic pulse parameters such as intensity and spec- 
tral width and to investigate dynamical effects. By com- 
paring projections and actual final state populations in 
the previous section, we have shown that a simple two- 
photon Franck-Condon principle does not correctly cap- 
ture the stabilization dynamics. A less simplified picture 
is obtained by taking the structure of the vibrational lev- 
els in the initial, intermediate and final electronic states 
fully into account, neglecting strong field effects and a 
dynamical interplay between pulse and spin-orbit cou- 
plings. Specifically, in the weak-field regime and for TL 
pulses, the probability of the resonant two-photon transi- 
tion is obtained by perturbation theory. It is determined 
by the effective two-photon transition moment, 



D{v") = 



E ■ 

j=2,3,4 



(12) 

where (vF m |d>') = £\ i# / ^\ n {R)*d j {R)^,{R)dR de- 
notes the electric transition dipole moment between the 
initial wave packet with components i in the 
(6) 1 S+, and (3) 3 n g manifold and the intermediate vibra- 
tional levels v' with components i' in the excited A 1 !!^ 
andfe 3 n g . Correspondingly, (v'\d\\v") is the electric tran- 
sition dipole moment between the intermediate levels v' 
of the A ll E+ and b 3 H g states and the final ground X 1 ^^ 
state vibrational level v" . The exponent accounts for the 
bandwidth of the laser pulse, and the detuning of the 
intermediate levels from the laser frequency, 



Auv = (J2 Pl /2 (5 P ) — w v i + u) v " — lol 



(13) 



The laser frequency, wl, is chosen such that the two- 
photon transition is resonant for the initial wave packet 
and the final level v", wi = (^p 1/2 (6p) ~ w m + w V ")/2. 
Here, w;„ denotes the binding energy of the initial wave 
packet, defined to be positive, and ^ 2 P 1/2 (n P ) is the exci- 
tation energy of the 2 -P1/2 atomic level. The spectral 
width is given in terms of the FWHM of the time profile 
of the pulse, r, er w = 2y2 ln(2)/r. 

The two-photon transition moments, defined by 
Eq. 12, are shown in Fig. 7 for all vibrational levels v" 
of the X 1 ^^ ground electronic state. The large peak 
around a binding energy of 2650 cm -1 indicates that sta- 
bilization of the excited state wave packet to levels with 
binding energies in this range is most efficient. The 
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FIG. 7: (Color online) Resonant two-photon transition mo- 
ments defined by Eq. (12) vs binding energy of the target 
X state levels v" for four different pulse durations or spec- 
tral widths, respectively. 



peak maximum in Fig. 7 corresponds to a transition to 
the vibrational level v" = 23, with a binding energy 
E v ti=2651 cm -1 . This suggests that the level v" = 23 
might be a good choice for the target in the state-to- 
state optimization of stabilization process below in Sec- 
tion V. A standard choice of 100 fs pulse duration for 
the TL pulse, corresponding to a spectral width of about 
150 cm" 1 FWHM is sufficient to address a broad distri- 
bution of target vibrational levels in the X 1 ^^ state. 
The integrated pulse energy is given by 



E = e cA / \e(t)\ dt, 
Jo 



(14) 



with e(t) the laser field, A = 7rr 2 the area which is cov- 
ered by the laser (r = 50 fim was assumed) , c the speed 
of light, and £o the dielectric constant. We use the in- 
tegrated pulse energy rather than the peak intensity of 
the pulse since, independently of the pulse duration, it 
quantifies the energy pumped into the molecule. 

The two-photon transition probability can be predicted 
from the effective two-photon transition moment, cf. 
Eq. (12) and Fig. 7, only in the weak-field regime when 
dynamic Stark shifts and other time-dependent effects 
do not play any role. The dependence of the two-photon 
transition probability on the pulse intensity and detun- 
ing, 



Awl = 2uj l -W2 P 



i(6p) 



(15) 



is illustrated in Fig. 8(a). The pulse duration is kept 
fixed at 100 fs FWHM. For weak and intermediate pulse 
intensities, with the integrated pulse energy correspond- 
ing to 8.3 nJ and 25.4 nJ, the final X 1 !]^ ground state 
population as a function of the pulse detuning reflects 
the shape of the effective two-photon transition moment, 
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FIG. 8: (Color online) a: Final population of the X 
state after a TL 100 fs pulse vs pulse detuning, taken with 
the respect to the 2 -Pi/2(6p) — s> 2 S(5s) atomic two-photon 
transition, for five different integrated pulse energies, b: 
Vibrational decomposition of the X 1 E+ component of the 
wave function after stabilization with TL 100 fs pulses with 
different pulse energies and detunings with the respect to 
2 -Pi/2(6p) — > 2 S(5s) atomic two-photon transition: □ - 
25.4nJ and 2650cm" 1 , ■ - 2.06^J and 2650cm" 1 , O - 
25.4nJ and 1240cm" 1 , • - 2.06 and 1240cm" 1 . 



Fig. 7. On the other hand, the final A" 1 E+ population 
decreases for the detuning corresponding to the maxi- 
mum of the two-photon transition moment and increases 
for smaller detunings when the integrated pulse energy 
is increased. This observation is rationalized in terms of 
the strong dynamic Stark shift by analyzing the vibra- 
tional distribution of the final X 1 ^ state population in 
Fig. 8(b) for a detuning, Awl=2650 cm" 1 , corresponding 
to the maximum of the two-photon transition moment in 
Fig. 7. When increasing the integrated pulse energy from 
25.4 nJ to 2 //J, i.e., from the intermediate to the strong 
field regime, two peaks are observed in the vibrational 
distribution rather than a single Gaussian around the 
binding energy of v" = 23, reflecting the bandwidth of 
the pulse. In the strong field regime, one peak of the 
vibrational distribution is still located around the bind- 
ing energy of the resonant level, while the second one 
is shifted by 800 cm" 1 to smaller binding energies. This 
is due to the positive differential Stark shift caused by 
the coupling to the intermediate states, which indeed in- 
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creases the energy separation between ground and ex- 
cited states by about 800 cm -1 . The dynamic Stark shift 
of the ground state vibrational level v" is estimated ac- 
cording to 



IKK 



UJ-, 



(16) 



where 



is the electric transition dipole mo- 



ment and u! V ' V " the transition frequency between levels 
v' and v", and ujl denotes the laser frequency. Fig- 
ure 8(b) compares the final state vibrational distribu- 
tion for two different dctunings of the stabilization pulse, 
Aojl = 2650 cm -1 corresponding to the peak of the two- 
photon transition probability for weak and intermediate 
fields (black dotted and red solid curves in Fig. 8(a)) 
and Aojl = 1240 cm -1 corresponding to the peak of 
the two-photon transition probability for strong fields 
(black dot-dashed curve in Fig. 8(a)). Inspection of 
the vibrational distributions in Fig. 8(b) reveals that for 
Aojl = 1240 cm -1 and 2 /iJ integrated pulse energy, a 
peak at binding energies larger than the detuning ap- 
pears. Also, this peak is caused by the differential Stark 
shift which this time is negative, decreasing the energy 
separation between ground and excited states by about 
600 cm -1 . In the weak and intermediate field regime, al- 
most no population is transferred for Alol = 1240 cm -1 
(red empty circles in Fig. 8(b)), confirming a strong field 
effect. 

The total population that is transferred by a TL 
pulse, with 100 fs FWHM and a detuning in the range 
of 2500 cm -1 to 2600 cm -1 , from the initial wave packet 
to the ground electronic state amounts to up to 

1.7% in the weak and intermediate field regime. For 
strong fields, up to 2.9% of the population can be trans- 
ferred for dctunings between 1000 cm" 1 and 1500 cm -1 
and pulse energies above 1 /^J. The subsequent analysis 
is restricted to pulses with detunings between 2500 cm -1 
and 2650 cm" 1 corresponding to the maximum of the ef- 
fective two-photon transition moment where the smallest 
pulse intensities should be required. 

As apparant from Fig. 8, the dynamic Stark shift is 
detrimental to efficient population transfer by the sta- 
bilization pulse. One option to increase the integrated 
pulse power while keeping the maximum field intensity, 
and thus the dynamic Stark shift, small is to consider a 
train of short TL pulses. A second option is given by 
chirping the pulse. The efficiency of the two-photon pop- 
ulation transfer to the electronic ground with the first 
option, a train of 100 fs TL pulses delayed relative to 
each other by 200 fs and with a sinusoidal envelope, is 
analyzed in Fig. 9. While increasing the pulse energy of 
a 100 fs pulse does not improve the population transfer to 
the electronic ground state beyond 1.7%, a train of pulses 
yields up to about 6% for pulse energies that are still in 
the nJ range. The population transfer with a train of 
short pulses can be interpreted as the cumulative result 
of many single transitions that accumulate amplitude in 
the X 1 ^^ ground state. Using a train of pulses instead 
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FIG. 9: (Color online) a: X l Y^ state population after a 100 fs 
TL pulse (solid line) and a train of 100 fs TL pulses (dashed 
and dotted lines) vs integrated pulse energy, b: state 
vibrational distribution after stabilization with a 100 fs TL 
pulse and 4ps train of 100 fs TL pulses. The pulse detuning 
is 2500cm" 1 taken with the respect to 2 A/2(6p) — > 2 S(5s) 
atomic two-photon transition. 



of a single pulse with the same bandwidth, one can pro- 
duce 3.5 times more ground state molecules. Fig. 9 thus 
confirms that the Stark shift is responsible for the com- 
paratively inefficient population transfer observed for TL 
pulses. 

From the coherent control of atomic transitions using 
strong fields, it is known that the influence of the dy- 
namic Stark shift can be compensated by chirping the 
pulse [14, 15]. We investigate in Fig. 10(a) showing the 
final ground state population vs chirp rate for increas- 
ing pulse energy whether this approach works also for 
molecular transitions. We use a positive chirp to correct 
the influence of the dynamic Stark shift since the dif- 
ferential Stark shift for stabilization to vibrational lev- 
els with binding energies close to 2650 cm -1 is positive. 
Chirping a weak-field pulse (black dot-dashed curve) de- 
teriorates the population transfer. When more energetic 
pulses are used, chirping increases the final X 1 ^^ state 
population from 1.5% for unchirped pulses to almost 9% 
for the best chirped pulses. In total we find that chirp- 
ing the pulse improves the stabilization process and en- 
hances the amount of ground state population by a factor 
of about six. Figure 10(b) showing the final state vibra- 
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FIG. 10: (Color online) a: Final X 1 ^ state population after 
linearly chirped pulses with a TL FWHM of 100 fs vs chirp 
rate for different pulse energies, b: X l Yi% state vibrational 
distribution after stabilization with a weak TL pulse, a TL 
and a chirped strong pulse. The pulse detuning is 2500 cm" 
taken with the respect to 2 A/2(6p) — > 2 S(5s) atomic two- 
photon transition. 



tional distribution confirms that the same mechanism as 
in the atomic case is at work [14, 15]: When increasing 
the pulse energy from 25.4 nJ to 300 nJ without chirping 
the pulse, a second peak shifted by 500 cm -1 appears. 
The energies of the levels of this second peak correspond 
exactly to the detuning corrected by the Stark shift. A 
linear chirp introduces a time-dependent instantaneous 
frequency of the pulse, uj(t) = lul + X^/2 ; with x the 
temporal chirp rate. When chosen correctly, the chirp 
compensates the phase that the molecules accumulate 
due the Stark shift and thus prevents the transition to 
shift out of resonance. This leads to the strong enhance- 
ment of the stabilization efficiency observed in Fig. 10. 

Our investigation of the stabilization dynamics under 
TL and linearly chirped pulses shows that simply re- 
placing a strong TL pulse by a train of pulses with the 
same total integrated pulse energy or linearly chirping 
the pulse can enhance the stabilization probability from 
1% up to 9%. The reason for the enhancement is given by 
the weaker Stark shifts for smaller peak intensities and 
compensation of the phase accumulated due to the Stark 
shift by a linear chirp. In the following section we will 
employ optimal control theory to calculate the optimum 
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FIG. 11: (Color online) a: Final X E^" state population af- 
ter optimized pulses of different pulse duration vs number of 
iteration steps, b: Final X £jj~ state population (left scale, 
black circles) and corresponding pulse energy (right scale, red 
squares) vs pulse duration. The values were collected after 
50 iterations, starting with the same weak-field guess pulse 
(E = 10 nJ) and using the same weight, A = 400, in all opti- 
mizations. 



detunings and pulses shapes. This will allow us to de- 
termine the maximum number of ground state molecules 
that can be produced for a given integrated pulse energy. 



V. STABILIZATION TO THE ELECTRONIC 
GROUND STATE WITH OPTIMALLY SHAPED 
PULSES 

We employ the optimization algorithm described in 
Sec. IIC to find those laser pulses that stabilize the 
initial wave packet most efficiently to the ^ X E+ elec- 
tronic state. The final X 1 E+ state population, shown in 
Fig. 11(a), converges smoothly to the maximal value that 
can be obtained with a given pulse duration, displayed 
in Fig. 11(b). The maximum stabilization probability 
for a pulse duration of 1 ps is 13.7%. Increasing the pulse 
duration, the stabilization probability reaches 26.4% for 
2 ps pulse and 67% for 10 ps pulses. The integrated pulse 
energies of the optimized pulses vary between 80 nJ and 
150 nJ. This is two to three times smaller than the inte- 
grated pulse energies for the trains of pulses and the lin- 
early chirped pulses discussed in Sec. IV. The guess pulse 
for the optimizations shown in Fig. 11 is a TL pulse with 
a pulse duration of 100 fs and integrated pulse energy 
of 10 nJ. For all the results presented here, the state-to- 
manifold-of-states functional, Eq. 9, was employed. The 
results obtained by using the state-to-state functional, 
Eq. 8, do not differ significantly. In particular, the same 
bounds on the maximum stabilization efficiencies are ob- 
served. However, for the state-to-state functional, the 
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optimizations were found to converge slower. This is eas- 
ily rationalized in terms of a single state being a more 
restrictive optimization target than a manifold. The in- 
tegrated energy of optimal pulses presented in Fig. 11(b) 
does not depend strongly on the optimal pulse duration. 
The slightly oscillatory behavior of the integrated pulse 
energy as a function of pulse duration is observed irre- 
spective of the shape and energy of the guess pulse and 
the weight A. 

The fact that the maximum population transfer to the 
X 1 E+ state is clearly less than 100% is due to the pulses 
being too short to drive the complete wave packet to 
the ground state [60]. When the pulse duration is much 
shorter than the time scale of the vibrational motion and 
spin-orbit oscillations, then only that part of the wave 
packet that shows a favorable overlap with the target 
state during the optimization window is transferred. For 
example, the pulse with 1 ps duration essentially reflects 
the overlap of the initial wave packet. By increasing the 
pulse duration, cf. Fig. 11(b), the stabilization efficiency 
increases monotonically. A stabilization probability of 
100% is expected once the pulse duration is longer than 
the revival time of the wave packet, estimated above to 
be between 20 ps and 30 ps. In fact, optimizations with 
pulse durations of 20 ps and 30 ps yield stabilization effi- 
ciencies of 93% and 99%, respectively, with pulse energies 
below 150 nJ. However, we restrict our analysis to pulse 
durations up to 10 ps since stretching a TL 100 fs pulse 
by more than a factor of 100 due to the pulse shaping is 
not realistic. 

Analyzing the time evolution of the population on each 
of the electronic states during an optimized pulse of 10 ps 
pulse duration, the molecules are found to first accumu- 
late in the intermediate and 6 3 n u states before be- 
ing dumped to the A x 5]+ electronic ground state. The 
example of an optimized pulse with pulse duration of 
8ps is presented in Fig. 12. Inspection of the optimized 
pulse during a short interval of 400 fs, Fig. 12(b), reveals 
that each peak of the pulse amplitude is correlated to 
a change of the temporal phase by n. The spectrum, 
Fig. 12(c), displays two pronounced peaks with maxima 
at ±400 cm -1 with respect to the central laser frequency. 
The disappearence of the central frequency during the 
optimization and appearance of two slightly detuned fre- 
quencies is somewhat surprising since the central fre- 
quency was chosen to maximize the effective two-photon 
transition moment. There are two possibilities to ratio- 
nalize this result of the optimization: Either the detuning 
shifted by the additional 400 cm -1 is better and should 
be chosen for the guess pulse, or the absorption of two 
photons with different energies is more optimal than that 
of two identical photons. The latter explains the ob- 
served optimal spectrum: Analysis of the electric transi- 
tion dipole moment between the initial wave packet and 
the intermediate vibrational levels, (^i n \d e \v') , reveals 
that it takes its maximal value for levels v' that are de- 
tuned from the frequency corresponding to the maximum 
of the effective two-photon transition moment, Eq. (12), 
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FIG. 12: (Color online) Temporal envelope (a) and spectrum 
(c) of an optimized pulse, b: Envelope (left scale, solid black 
line) and temporal phase (right scale, red dashed line) of the 
optimized pulse in a short time interval. 



by 460 cm -1 . The electric transition dipole moment be- 
tween the intermediate and the ground state vibrational 
level, (v'|d e |u"), has its maximum for a transition fre- 
quency smaller by 250 cm" 1 than the laser frequency cor- 
responding to the maximum of the effective two-photon 
transition moment. Note that the transition moments 
for the absorption of the first photon, (^ !rl |de|w'), are 
about 7 times smaller than those for the absorption of 
the second photon, (v'\d e \v"). The effective two-photon 
transition moment that is obtained as a compromise of 
the two one-photon transition moments, cf. Eq. (12), has 
its maximum at a detuning of 2650 cm -1 , cf. Fig. 12. 
Allowing for two photons of different energy in the cal- 
culation of the effective two- photon transition moment, 
we still find a peak for a ground state binding energy of 
2650 cm -1 , which is at best 40% higher when energies of 
photons are detuned by ±390 cm -1 . It corresponds to 
the transition frequency from the initial wave packet to 
the intermediate state being 390 cm -1 larger and that be- 
tween the intermediate state and the ground electronic 
state being 390 cm -1 smaller than the frequency for a 
transition with two identical photons. Equipped with 
this information, we can construct a guess pulse that is 
the sum of two TL pulses with their central frequencies 
separated by 780 cm -1 . In this case, half of the inte- 
grated guess pulse energy is sufficient to reach the same 
initial stabilization probability, reflecting the stronger 
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effective two-photon transition moment. Optimization 
with such a guess puise converges faster and the final in- 
tegrated pulse energy of the optimized pulse is smaller 
(data not shown) but the bound for the stabilization ef- 
ficiency, Fig. 11, remains in place. The spectrum of the 
optimized pulse with two peaks separated by 780 cm -1 
is very similar to the one shown in Fig. 12, irrespective 
of the guess pulse central frequency. 

Optimization of the stabilization pulse reveals that the 
upper bound of the stabilization efficiency, found to be 
9% for linearly chirped pulses, can be increased up to 
67% when a TL pulse of 100 fs pulse duration is shaped 
and stretched to 10 ps. At the same time, the integrated 
pulse energies of the optimized pulses are below 150 nJ, 
two to three times less than those found for the best 
linearly chirped pulses or trains of TL pulses in Sec. IV. 

VI. SUMMARY AND CONCLUSIONS 

Based on state of the art ab initio calculations, we 
have studied the optical production of Rb2 molecules in 
their electronic ground state using multi-photon transi- 
tions that are driven by short laser pulses. Our model 
includes not only accurate potential energy surfaces but 
also spin-orbit couplings, Stark shifts and transition ma- 
trix elements. We have employed a non-perturbative 
treatment of the light-matter interaction which is crucial 
to capture the strong-field effects that often accompany 
multi-photon transitions. 

Our proposal for the optical production of molecules 
using shaped femtosecond laser pulses that drive multi- 
photon transitions consists, in its first step, of non- 
resonant three-photon photoassociation of atom pairs 
colliding in their triplet state. A three-photon transition 
allows to access electronic states that vary as 1/R 3 at 
long range, providing comparatively large frec-to-bound 
transition matrix elements [22]. Strong spin-orbit inter- 
action allows for triplet-to-singlct conversion. The sta- 
bilization pulse, time-delayed with respect to the pho- 
toassociation pulse, transfers the photoassociated wave 
packet to the electronic ground state in a resonant two- 
photon transition proceeding via the ^4 1 E+-6 3 n u man- 
ifold. It benefits from the intricate excited state wave 
packet dynamics resulting from coupled vibrational dy- 
namics and singlet-triplet oscillations due to the spin- 
orbit interaction. 

We have studied the transfer of the excited state wave 
packet to the electronic ground state using transform- 
limited, linearly chirped and optimally shaped laser 
pulses. Surprisingly, we find the integrated pulse en- 
ergy of the optimally shaped laser pulses to be signif- 
icantly lower than that of the best transform-limited 
and linearly chirped pulses, while yielding a much bet- 
ter stabilization efficiency. This is due to the fact that 
the transform-limited and linearly chirped laser pulses 



were chosen based on an effective two-photon transition 
matrix clement assuming equal transition frequencies of 
both photons. Optimization reveals that a two-photon 
transition with two slightly different transition frequen- 
cies allows to employ two single-photon transitions with 
significantly larger transition matrix elements. Overall, 
the stabilization efficiency is limited by somewhat less 
than 70% for transform-limited 100 fs laser pulses that 
are stretched to 10 ps. More than 90% transfer efficiency 
becomes possible by stretching the pulse to 20 ps. 

Linearly chirped pulses were found to perform much 
better, by almost an order of magnitude, than transform- 
limited pulses. This is due to large Stark shifts which 
drive the transition off resonance for transform limited 
pulses. We have confirmed that a strong-field control 
scheme known for atomic transitions, with a linear chirp 
compensating the phase accumulated due to the Stark 
shift [14, 15], can also be successfully employed for molec- 
ular transitions. In this case, the chirp rate cannot be 
calculated analytically but needs to be determined nu- 
merically. 

Multi-photon transitions allow for utilizing the full 
bandwidth of a femtosecond pulse for photoassociation 
and stabilization. They also provide more flexibility 
with respect to the transition frequencies that can be ad- 
dressed. The basic features of femtosecond laser pulses - 
their broad bandwidth and pulse shaping capabilities - 
can then be used to full advantage. Our results thus pave 
the way towards the coherent control of photoassociation 
at very low temperature. They are not restricted to ru- 
bidium but rather are applicable to any molecule that 
shows the main features of the dynamics studied here - 
an ion-pair potential that is coupled to covalent states 
and strong spin-orbit interaction. 

Future work will consider optimization of the pho- 
toassociation pulse. This represents on non-trivial con- 
trol problem since the initial state of photoassociation is 
the thermally populated continuum of scattering states. 
While this thermal ensemble is inherently incoherent, 
quantum effects are already perceptible due to the low 
temperature. It will be interesting to see whether an 
optimally shaped pulse can make use of the enhanced 
quantum purity and presence of correlations at low tem- 
perature. 



Acknowledgments 

We would like to thank Ronnie Kosloff for many 
fruitful discussions. Financial support from the Polish 
Ministry of Science and Higher Education through the 
project N N204 215539 is gratefully acknowledged. MT 
is supported by the project operated within the Founda- 
tion for Polish Science MPD Programme co-financed by 
the EU European Regional Development Fund. 



14 



K. M. Jones, E. Tiesinga, P. D. Lett, and P. S. Julienne, 
Rev. Mod. Phys. 78, 483 (2006). [28 
S. A. Rice and M. Zhao, Optical control of molecular 
dynamics (John Wiley & Sons, 2000). [29 
P. Brumer and M. Shapiro, Principles and Applications 
of the Quantum Control of Molecular Processes (Wiley 
Interscience, 2003). [30 
C. P. Koch and M. Shapiro, Chem. Rev. (2012). 

C. P. Koch, R. Kosloff, and F. Masnou-Seeuws, Phys. 

Rev. A 73, 043409 (2006). [31 
W. Salzmann, U. Poschinger, R. Wester, M. Wei- 
demiiller, A. Merli, S. M. Weber, F. Sauer, M. Plewicki, [32 

F. Weise, A. Mirabal Esparza, et al., Phys. Rev. A 73, 
023414 (2006). [33 

B. L. Brown, A. J. Dicks, and I. A. Walmsley, Phys. Rev. 
Lett. 96, 173002 (2006). [34 
W. Salzmann, T. Mullins, J. Eng, M. Albert, R. Wester, 

M. Weidemiiller, A. Merli, S. M. Weber, F. Sauer, 

M. Plewicki, et al., Phys. Rev. Lett. 100, 233003 (2008). [35 

T. Mullins, W. Salzmann, S. Gotz, M. Albert, J. Eng, 

R. Wester, M. Weidemiiller, F. Weise, A. Merli, S. M. [36 

Weber, et al., Phys. Rev. A 80, 063416 (2009). 

A. Merli, F. Eimer, F. Weise, A. Lindinger, W. Salz- [37 

mann, T. Mullins, S. Gotz, R. Wester, M. Weidemiiller, 

R. Aganoglu, et al., Phys. Rev. A 80, 063417 (2009). [38 

D. J. McCabe, D. G. England, H. E. L. Martay, M. E. 
Friedman, J. Petrovic, E. Dimova, B. Chatel, and I. A. [39 
Walmsley, Phys. Rev. A 80, 033404 (2009). 

D. Meshulach and Y. Silberberg, Nature 396, 239 (1998). [40 
D. Meshulach and Y. Silberberg, Phys. Rev. A 60, 1287 [41 
(1999). [42 

C. Trallero-Herrero, D. Cardoza, T. C. Weinacht, and [43 
J. L. Cohen, Phys. Rev. A 71, 013423 (2005). 

C. Trallero-Herrero, J. L. Cohen, and T. Weinacht, Phys. [44 
Rev. Lett. 96, 063603 (2006). 

C. Trallero-Herrero and T. C. Weinacht, Phys. Rev. A [45 
75, 063401 (2007). 

A. Gandman, L. Chuntonov, L. Rybak, and Z. Amitay, [46 
Phys. Rev. A 75, 031401 (2007). 

A. Gandman, L. Chuntonov, L. Rybak, and Z. Amitay, 
Phys. Rev. A 76, 053419 (2007). [47 
L. Chuntonov, L. Rybak, A. Gandman, and Z. Amitay, 
Phys. Rev. A 77, 021403 (2008). 

L. Chuntonov, L. Rybak, A. Gandman, and Z. Amitay, [ 
J. Phys. B 41, 035504 (llpp) (2008). 

Z. Amitay, A. Gandman, L. Chuntonov, and L. Rybak, [49 
Phys. Rev. Lett. 100, 193002 (2008). 
C. P. Koch, M. Ndong, and R. Kosloff, Faraday Disc. [50 
142, 389 (2009). 

L. Rybak, S. Amaran, L. Levin, M. Tomza, R. Moszynski, 

R. Kosloff, C. P. Koch, and Z. Amitay, Phys. Rev. Lett. 

107, 273001 (2011). [51 

L. Rybak, Z. Amitay, S. Amaran, R. Kosloff, M. Tomza, 

R. Moszynski, and C. P. Koch, Faraday Disc. 153, 383 

(2011). 

R. Kosloff, S. Rice, P. Gaspard, S. Tersigni, and D. Tan- 

nor, Chem. Phys. 139, 201 (1989). [52 

T. Ban, H. Skenderovic, R. Beuc, I. K. Bronic, 

S. Rousseau, A. Allouche, M. Aubert-Frecon, and [53 

G. Pichler, Chem. Phys. Lett. 345, 423 (2001). 

T. Ban, R. Beuc, H. Skenderovi, and G. Pichler, Euro- [54 



phys. Lett. 66, 485 (2004). 

C. Amiot, O. Dulieu, and J. Verges, Phys. Rev. Lett. 83, 
2316 (1999). 

C. M. Dion, C. Drag, O. Dulieu, B. Laburthe Tolra, 
F. Masnou-Seeuws, and P. Pillet, Phys. Rev. Lett. 86, 
2253 (2001). 

H. K. Pechkis, D. Wang, Y. Huang, E. E. Eyler, P. L. 
Gould, W. C. Stwalley, and C. P. Koch, Phys. Rev. A 
76, 022504 (2007). 

A. Fioretti, O. Dulieu, and C. Gabbanini, J. Phys. B 40, 
3283 (2007). 

S. Ghosal, R. J. Doyle, C. P. Koch, and J. M. Hutson, 
New J. Phys. 11, 055011 (2009). 

B. E. Londono, J. E. Mahecha, E. Luc-Koenig, and 
A. Crubellier, Phys. Rev. A 80, 032511 (2009). 

M. Tomza, F. Pawlowski, M. Jeziorska, C. P. Koch, 
and R. Moszynski, Phys. Chem. Chem. Phys. 13, 18893 
(2011). 

J. M. Sage, S. Sainis, T. Bergeman, and D. DeMille, 
Phys. Rev. Lett. 94, 203001 (2005). 

M. Ndong and C. P. Koch, Phys. Rev. A 82, 043437 
(2010). 

V. Kokoouline, O. Dulieu, R. Kosloff, and F. Masnou- 
Seeuws, J. Chem. Phys. 110, 9865 (1999). 
K. Willner, O. Dulieu, and F. Masnou-Seeuws, J. Chem. 
Phys. 120, 548 (2004). 

S. Kallush and R. Kosloff, Chem. Phys. Lett. 433, 221 

(2006) . 

R. Kosloff, Annu. Rev. Phys. Chem. 45, 145 (1994). 
M. Tomza, M. Musial, and R. Moszynski, in preparation. 
M. Musial, J. Chem. Phys. 136 (2012). 

D. I. Lyakh, M. Musial, V. F. Lotrich, and R. J. Bartlett, 
Chem. Rev. 112, 182 (2012). 

R. J. Bartlett and M. Musial, Rev. Mod. Phys. 79, 291 

(2007) . 

T. G. A. Heijmen, R. Moszynski, P. E. S. Wormer, and 

A. van der Avoird, Mol. Phys. 89, 81 (1996). 

A. Marte, T. Volz, J. Schuster, S. Diirr, G. Rempe, 

E. G. M. van Kempen, and B. J. Verhaar, Phys. Rev. 
Lett. 89, 283202 (2002). 

R. Gutterres, C. Amiot, A. Fioretti, G. Gabbanini, 
M. Mazzoni, and O. Dulieu, Phys. Rev. A 66, 024502 
(2002). 

M. Marinescu and A. Dalgarno, Phys. Rev. A 52, 311 
(1995). 

I. S. Lim, P. Schwerdtfeger, B. Metz, and H. Stoll, J. 
Chem. Phys. 122, 104103 (2005). 

J. F. Stanton, J. Gauss, S. A. Perera, J. D. Watts, 
A. Yau, M. Nooijen, N. Oliphant, P. Szalay, W. Laud- 
erdale, S. Gwaltney, et al., ACES II is a program product 
of the Quantum Theory Project, University of Florida. 
H.-J. Werner, P. J. Knowles, F. R. M. R. Lindh, 
M. Schiitz, P. Celani, T. Korona, G. Rauhut, R. D. 
Amos, A. Bernhardsson, A. Berning, et al., MOLPRO, 
version 2006.1, a package of ab initio programs (2006), 
see http://www.molpro.net. 

S. J. Park, A. W. Suh, Y. S. Lee, and G.-H. Jeung, J. 
Mol. Spec. 207, 129 (2001). 

D. Edvardsson, S. Lunell, and C. M. Marian, Mol. Phys. 
101, 2381 (2003). 

A.-R. Allouche and M. Aubert-Frecon, J. Chem. Phys. 



15 



136, 114302 (2012). 
[55] J. P. Palao and R. Kosloff, Phys. Rev. A 68, 062308 
(2003). 

[56] S. E. Sklarz and D. J. Tannor, Phys. Rev. A 66, 053619 
(2002). 

[57] D. M. Reich, M. Ndong, and C. P. Koch, J. Chem. Phys. 
136, 104103 (2012). 



[58] W. Skomorowski, R. Moszynski, and C. P. Koch, Phys. 

Rev. A 85 (2012). 
[59] I. S. Averbukh and N. F. Perelman, Phys. Lett. A 139, 

449 (1989). 

[60] C. P. Koch, J. P. Palao, R. Kosloff, and F. Masnou- 
Seeuws, Phys. Rev. A 70, 013402 (2004). 



